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SUMMARY 


The  results  of  experiments  conducted  to  d  termine  the  dynamic  stability 
characteristics  of  a  dynamically  similar  model  of  a  quad-duct  v/STOL  air¬ 
craft  are  reported  in  References  1  and  2.  Portions  of  these  data  that 
pertain  to  the  longitudinal  dynamics  of  the  vehicle  at  five  duct  incidences 
were  analyzed  to  determine  the  stability  derivatives  of  the  vehicle.  The 
analysis  and  the  resulting  stability  derivatives  are  presented  and  dis¬ 
cussed  in  this  report. 

The  measured  time  histories  indicated  that  the  data  could  be  analyzed  on 
the  basis  of  linearized  small  perturbation  equations.  Root  locus  tech¬ 
niques  were  used  to  analyze  the  data. 

The  full-scale  derivatives  determined  from  the  analysis  that  correspond  to 
a  ve  icle  very  similar  to  the  Bell  X-22A  are  presented. 

The  transient  motions  of  the  model  were  unstable  at  all  duct  incidences 
except  50°,  the  lowest  incidence  investigated. 
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INTRODUCTION 


References  1  and  2  present  data  from  a  series  of  experiments  conducted  to 
measure  the  longitudinal  transient  response  characteristics  of  a  dynami¬ 
cally  similar  quad-duct  V/STOL  model  using  the  Princeton  Dynamic  Model 
Track.  This  program  is  part  of  a  continuing  effort  using  the  Princeton 
Dynamic  Model  Track  to  provide  data  on  the  dynamic  stability  of  v/STOL 
aircraft  at  low  speeds.  The  model  is  similar  in  configuration  to  the  Bell 
X-22A,  with  a  scale  factor  of  0.145.  The  model  is  shown  mounted  on  the 
test  apparatus  in  Figure  1.  A  close-up  view  of  the  model  is  shown  in 
Figure  2.  A  general  arrangement  drawing  of  the  model  is  given  in  Figure  3. 
Differences  between  the  Bell  X-22A  and  this  research  model  are  described 
in  the  section  entitled  Description  of  Apparatus  and  Experiments. 

This  report  analyzes  the  data  from  References  1  and  2  to  determine  the 
longitudinal  stability  derivatives  of  the  vehicle  in  low- speed/high-duct  - 
incidence  flight  conditions,  including  hovering  flight.  The  test  program 
is  described  in  detail  in  References  1  and  2;  thus,  only  a  brief  discussion 
of  the  program  is  included  in  this  report.  The  test  conditions  analyzed 
are  listed  in  Table  I. 

Transient  measurements  in  the  flight  conditions  of  interest  indicated  that 
the  dynamic  motions  of  the  model  could  be  described  by  conventional, 
linearized  small-perturbation  equations.  The  following  analysis  is  based 
on  that  assumption.  Discussion  of  the  analysis  of  the  data  is  phrased  in 
terms  of  model  parameters.  The  values  of  the  stability  derivatives  of 
the  full-scale  aircraft  are  discussed  in  the  section  entitled  Stability 
Derivatives  of  the  Full-Scale  Aircraft.  The  scaling  laws  used  to  design 
the  model  result  in  the  conversion  factors  given  in  Table  II,  which  are 
used  for  interpretation  of  model  data  in  full-scale  terms. 

The  model  was  found  to  be  dynamically  unstable  at  all  but  the  lowest  duct 
incidence  examined.  Data  were  taken  on  the  response  of  the  basic  model, 
as  well  as  with  varying  amounts  of  rate  feedback  (fore  and  aft  differential 
propeller  blade  angle  proportional  to  pitching  rate),  to  assist  in  the 
analysis. 

It  will  be  noted  that  in  the  flight  conditions  at  duct  incidences  of  60° 
and  50° ,  the  eg  of  the  model  was  displaced  from  the  pivot  axis  of  the 
model.  Balance  weights  were  added  to  the  model  at  these  trim  conditions 
to  reduce  the  amount  of  differential  propeller  blade  angle  required  for 
pitching  moment  equilibrium. 

For  analysis  purposes,  it  is  possible  to  locate  the  origin  of  the  axis 
system  either  at  the  eg  of  the  model  or  at  the  pivot  axis.  It  was  con¬ 
sidered  to  be  more  convenient  to  locate  the  origin  of  the  axis  system  at 
the  pivot  axis  and  to  add  terms  to  account  for  the  displacement  of  the  eg 
from  this  location.  The  stability  derivatives  are  presented  about  the 
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pivot  axis,  which  corresponds  to  the  eg  position  of  the  Bell  X-22A  given 
in  Reference  3  (WL  139,  FS  312).  All  linear  velocities  of  the  model  are 
the  linear  velocities  of  the  pivot  axis. 

Various  methods  of  analyzing  the  data  are  possible.  A  convenient  approach, 
when  limited-degree-of-freedom  data  are  available,  employs  root  locus 
techniques  as  described  in  the  section  entitled  Experimental  Results  and 
Analysis  of  Data. 

While  only  the  frequency  and  damping  characteristics  of  the  motions  are 
presented  and  used  in  the  analysis,  various  other  properties  of  the  time 
histories  may  be  used  to  provide  additional  information.  Other  measured 
data,  such  as  the  phase  angle  between  the  pitch  angle  and  the  horizontal 
velocity  perturbation,  were  used  to  check  the  results  obtained  from  the 
method  of  analysis  described  in  the  text  that  follows. 
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DESCRIPTION  OF  APPARATUS  AND  EXPERIMENTS 


TEST  F'CILITY 

The  Princeton  Dynamic  Model  Track  is  designed  expressly  for  the  study  of 
the  dynamic  motions  of  helicopter  and  V/STOL  moaels  at  equivalent  flight 
speeds  of  up  to  60  knots  (for  a  one-tenth  scale  model).  Basic  components 
of  the  facility  include  a  servo-driven  carriage  riding  on  a  track  750  feet 
long,  located  in  a  building  with  a  cross  section  of  30  by  30  feet;  the 
carriage  has  an  acceleration  potential  of  0.6g  and  a  maximum  speed  of  40 
feet  per  second.  A  detailed  description  of  uhe  facility  and  the  testing 
tecnniques  employed  may  be  found  in  Reference  4. 

A  model  can  be  attached  to  the  carriage  by  one  of  several  bocms.  The 
mount  used  to  conduct  longitudinal  investigations  is  shown  in  Figure  1. 
This  mount  permits  relative  displacements  of  the  model  with  respect  to 
the  carriage  in  horizontal  and  vertical  directions.  The  model  is  sup¬ 
ported  on  a  three-axis  gimbal  system  that  allows  selection  of  any  or  all 
of  the  three  angular  degrees  of  freedom.  Horizontal  relative  motion  of 
the  model  with  respect  to  the  carriage  is  sensed  and  used  to  command  tne 
carriage  to  follow  the  model  in  a  closed-loop  fashion  Similarly,  vertica 
displacement  of  the  model  with  respect  to  the  carriage  commands  the  boom 
to  move  vertically.  This  servo  operation  of  the  carriage  allows  the  model 
to  fly  "free",  with  no  restraints  on  the  dynamic  motions  being  investi¬ 
gated.  This  method  of  testing  may  be  considered  to  be  similar  to  dynamic 
flight  testing,  but  considerably  more  control  over  the  experiment  is 
possible . 

MODEL 


A  photograph  of  the  model  is  shown  in  Figure  2.  A  general  arrangement 
drawing  is  presented  in  Figure  3-  Figure  4  shows  the  pivot  axis,  eg 
location,  and  reference  locations  of  the  model.  The  model's  pertinent 
dimensions  and  inertia  characteristics  are  listed  in  Table  III.  Tne  model 
was  designed  as  a  general  research  model  for  invest,! gation  of  the  dynamic 
stability  characteristics  of  various  quad  configuration  V/STOL  aircraft  as 
described  in  Reference  5. 

This  dynamic  model  is  powered  by  a  200-volt,  NOO-cycae ,  3-phase  electric 
motor.  The  motor  drives  the  four  ducted  propellers  through  a  central 
transmission  and  various  right-angle  gearboxes.  The  aerodynamic  shape  of 
the  model  is  obtained  through  the  use  of  a  Fiberglas  skin  with  Styrofoam 
stiffeners.  The  propeller  blades  are  made  with  a  plastic  foam  core  and 
Fiberglas  skin.  The  geometric  characteristics  of  the  propeller  are  mown 
in  Figure  5?  and  the  geometric  characteristics  of  the  duct  and  elevon  are 
shown  in  Figures  6  and  7.  The  duct  shape  is  identical  to  that  of  the 
Bell  X-22A  aircraft. 
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Nodal  control  positions  are  set  from  a  control  console  on  the  carriage. 

The  blade  pitch  angles  of  each  of  the  four  propellers  are  electrically 
controllable.  Also,  the  deflection  angles  of  the  elevons  are  electrically 
oontrollable.  All  of  these  control  systems  are  closed-loop  position 
controls  and  are  used  as  such  In  the  portions  of  the  experiments  involving 
feedback  to  alter  the  transient  motions  of  the  model.  The  dynamic  charac¬ 
teristics  of  these  feedback  loops  are  such  that  the  time  response  of  the 
control  is  negligible  in  the  frequency  range  of  interest.  Although  the 
control  servo  loops  are  nonlinear,  using  polarized  relays  for  power  ampli¬ 
fication,  they  can  be  characterized  as  having  a  closed-loop  natural 
frequency  of  approximately  10  cycles  per  second,  with  a  damping  ratio  of 
approximately  seven-tenths.  The  servo  gear  ratios  were  selected  so  that 
the  rate  limits  arising  from  the  rpm  limitations  of  the  control  drive 
motors  were  equal  to,  or  greater  than,  scaled  rate  limits  determined  from 
full-scale  ’Jell  X-22A  values,  as  given  in  Reference  3. 

This  research  model  differs  from  the  Bell  X-22A  in  the  following 
particulars: 

1.  The  elevon  on  the  model  differs  from  that  on  the  full-scale  air¬ 
craft.  The  model  elevon  has  no  movable  surface  forward  of  the 
hinge  line,  and  its  hinge  line  is  located  below  the  trailing  fdge 
of  the  duct,  as  shown  in  Figure  7.  While  these  differences  would 
affect  the  control  effectiveness  and  the  control  loads,  they, 
would  not  be  expected  to  have  any  significant  effect  on  the 
dynamic  motions. 

2.  The  duct  rotation  point  ia  at  a  different  location  on  the  model 
(8k  percent  c)  than  on  the  full-scale  aircraft  (55  percent  c). 

With  the  ducts  at  90°  incidence,  the  propeller  hubs  are  in  the 
same  relative  position  on  the  model  as  they  are  on  the  full-scale 
aircraft.  The  eg  of  the  model  is  higher  (1.2  percent  c)  on  the 
modal,  with  respect  to  the  propeller  hubs,  than  it  is  on  the 
full-scale  aircraft. 

3.  For  the  tests  at  90°,  80°,  and  70°  incidence,  the  vertical  tail 
on  the  model  is  smaller  than  the  one  on  the  full-scale  vehicle, 
as  shown  in  Figure  3>  The  scaled  vertical  tail,  also  shown  in 
Figure  3,  was  used  for  the  tests  at  duct  incidences  of  60°  and 
50®.  This  difference  in  vertical  tail  area  would  not  have  a 
significant  effect  on  the  longitudinal  dynamic  stability  charac¬ 
teristics. 

This  model  was  planned  as  a  general  research  model;  numerous  other  quad 
configuration  layouts  can  be  simulated  through  the  use  of  interchangeable 
parts  aa  described  in  Reference  5.  No  attempt  was  made  in  the  design 
stage  to  simulate  the  X-22A  precisely.  However,  the  modifications  de¬ 
scribed  above  will  not  result  in  appreciable  differences  in  the  model 
dynamic  stability  characteristics. 


k 


« 


EXPERIMENTS 


The  experiments  that  were  conducted  to  determine  the  stability  character¬ 
istics  of  the  model  consisted  of  transient  response  measurements  in 
various  longitudinal  degrees  of  freedom.  The  data  are  presented  in 
References  1  and  2,  and  the  test  conditions  are  given  in  Table  I.  Measure¬ 
ments  were  also  made  with  various  levels  of  rate  feedback.  Since  the 
model  was  unstable,  or  at  best  neutrally  stable,  in  the  majority  of  flight 
conditions  investigated,  no  predetermined  control  inputs  were  used  to 
excite  the  motion  of  the  model. 
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EXPERIMENTAL  RESULTS  AND  ANALYSIS  OF  DATA 


The  transient  response  data  presented  in  References  1  and  2  were  analyzed 
to  determine  the  period  and  damping  of  the  longitudinal  modes  of  motion  of 
the  model  at  five  duct  incidences;  90°,  80°,  70°,  60°,  and  50°.  The  trim 
condition  was  level  flight,  and  the  fuselage  attitude  was  set  equal  to 
zero  in  trim.  Experimentally  determined  model  trim  conditions  are  shown 
in  Figure  8.  The  data  and  the  results  of  the  analysis  are  discussed  in 
terms  of  model  parameters  in  this  section.  A  detailed  numerical  example 
of  the  analysis  procedure  at  a  duct  incidence  of  70°  is  presented  in 
Appendix  II. 

First,  it  is  desirable  to  nuike  p.  few  remarks  regarding  the  analysis  of 
transient  response  data  for  higher  than  second-order  systems.  If  a  system 
has  an  unstable  mode  present  in  its  transient  response,  it  is  difficult  to 
measure  (from  the  time  history  of  the  motion)  the  characteristics  of  any 
mode  but  the  unstable  one,  which  will  dominate  the  motion  irrespective  of 
the  nature  of  the  disturbance.  Therefore,  certain  practical  limitations 
are  placed  on  the  determination  of  all  of  the  stability  derivatives  of  an 
unstable  aircraft,  since  the  characteristics  of  all  the  modes  cannot  be 
accurately  determined.  This  limitation  in  analyzing  transient  response 
data  of  multiple-degree-of-freedom  systems  can  be  surmounted  by  the  use  of 
limited-degree -of-freedom  tests,  as  may  be  seen  from  the  following  dis¬ 
cussion. 

The  measured  characteristics  of  the  transient  response  consist  of  the 
frequency  and  damping  characteristics  of  the  dominant  unstable  mode.  When 
a  component  of  the  transient  motion  is  purely  divergent,  it  has  been  found 
difficult  to  measure  the  value  of  the  positive  real  root  corresponding  to 
the  divergence.  Therefore,  transient  response  data  which  involve  a  di¬ 
vergent  mode  have  not  been  analyzed  quantitatively.  Divergent  motions 
were  eliminated  in  the  single-degree-of-freedom  tests  (0  only)  through  the 
use  of  a  mechanical  spring.  The  test  conditions  in  which  a  pure  divergence 
was  present  were  two-degree -of- freedom  motions  (0-wf)  at  duct  incidences 
80°,  70°,  and  60°. 

Four  different  combinations  of  degrees  of  freedom  were  measured  at  all  of 
the  duct  incidences  except  hover:  one  single-degree-of-freedom  motion, 

0,  (Uf  =  0,  and  w^.  =  0);  two  two-degree-of-freedom  motions,  G-w^., 

(Uj.  =  0),  and  0-Uf-,  (w^,  =  0'S;  and  the  complete  longitudinal  three-degree- 

of-freedom  motion  (0-Uf-w^>'' .  The  development  of  the  equations  of  motion 

that  are  assumed  to  apply  to  the  data  is  presented  in  Appendix  I.  Note 
that  the  following  discussion  is  phrased  in  terms  of  space-fixed  degrees 
of  freedom  in  accordance  with  the  manner  in  which  the  tests  were  conducted. 
The  space-fixed  axis  system  is  shown  in  Figure  9>  and  the  transformation 
from  a  stability  axis  system  to  a  space-fixed  axis  system  is  discussed  in 
Appendix  I. 
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Measured  transient  characteristics  are  presented  at  the  various  duct  inci¬ 
dences  on  the  complex  plane  in  Figure  10.  Frequency  and  damping  charac¬ 
teristics  of  the  dominant  mode  are  presented  as  a  function  of  rate  feed¬ 
back  gain  and  degrees  of  freedom  of  the  test.  The  characteristic  roots  of 
the  single-degree-of-freedom  motion  are  presented  as  calculated  with  the 
mechanical  spring  removed  (see  Appendix  II ). 

The  following  general  approach  was  used  to  determine  the  stability  deriva¬ 
tives.  To  simplify  the  discussion,  known  additional  terms  due  to  the 
displacement  of  the  eg  of  the  model  from  the  pivot  axis  and  due  to  the 
model  mounting  linkage  mass  are  not  included  in  the  discussion  that 
follows  (see  Appendix  I). 

ONE  DEGREE  OF  FREEDOM 


First,  the  one-degree-of-freedom  results  were  analyzed.  At  all  duct  inci¬ 
dences  except  i^  =  50°,  a  mechanical  spring  was  added  about  the  pitch  axis 

of  the  model  to  provide  a  restoring  moment  proportional  to  model  attitude 
such  that  the  angular  motion  of  the  model  would  be  oscillatory,  and  thus 
could  be  analyzed  more  accurately  in  the  flight  conditions  where  the 
angle-of-attack  stability  of  the  model  was  positive  (M^Uq  >  0).  Positive 

angle -of- attack  stability  will  result  in  a  divergent  motion  of  the  model 
in  one  degree  of  freedom  if  no  restoring  spring  is  provided.  The  spring 
constant  of  the  mechanical  spring  was  selected  so  that  the  frequency  of 
motion  in  one  degree  of  freedom  was  similar  to  the  frequency  of  the  free 
motion  of  the  model  in  the  multiple-degree-of-freedom  tests.  At  the 
highest  speed  (lowest  duct  incidence)  tested,  the  angle-of-attack  stability 
was  negative  (M^U^  <  0)  and  therefore  no  spring  was  used. 

As  shown  in  Appendix  I,  the  equation  of  motion  that  applies  to  the  single- 
degree-of-freedom  tests  is  equation  (l8),  with  zcg  =  0. 

«  -  (»e  +  W>f)  ®  +(t~  '  VoJ e  ‘  0  (1> 


k0  is  the  spring  constant  of  the  mechanical  spring  and  is  determined  by 
em 

calibration  prior  to  the  experiment.  The  two  terms  in  equation  (l),  the 

coefficients  of  0  and  0,  are  determined  by  the  frequency  (u>N)  and  damping 

(ct)  measured  from  the  transient  response  data  (listed  in  Table  I).  De¬ 
rivatives,  or  combinations  of  derivatives,  are  found  from  the  relation¬ 
ships 
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+  Me  +  Viw{Jof  a  2p 


(2) 


Values  of  the  derivatives  calculated  from  the  measured  characteristics  of 
the  transients  and  the  expressions  of  equation  (2)  are  listed  in  Table  IV. 

TWO  DEGREES  OF  FREEDOM 

The  stability  derivatives  obtained  from  the  single-degree-of-freedom  runs 
are  now  used  in  conjunction  with  the  data  from  the  two-degree-of-freedom 
experiments  involving  pitch  angle  and  horizontal  velocity  to  find  other 
stability  derivatives.  The  equations  which  describe  the  motion  in  this 
case  are  equations  (16),  given  in  Appendix  I,  with  m/m^  equal  to  one  and 

zcg  equal  to  zero. 


uf  -  XyUf  +  (g  -  VJQ  )  0  =  0 

-  Vf  +  e  -  (m$  +  M*u0f )  e  -  »y;0f  e  =  o  (3) 


While  three  characteristic  roots  are  necessary  to  delineate  this  dynamic 
condition,  only  the  two  which  determine  the  oscillatory  mode  character¬ 
istics  can  be  evaluated,  as  discussed  previously.  There  are  three 
additional  derivatives  present  in  equations  (3):  Xy,My,  and  As  may 

be  seen  from  the  discussion  that  follows,  it  is  possible  to  calculate  only 
the  term  Xy  in  combination  with  My.  Analysis  of  isolated  duct  data  indi¬ 
cates  that  the  term  X^Uq  is  negligible  compared  to  g. 

A  convenient  way  of  calculating  the  two  unknown  derivatives  Xu  and  My  is 
to  place  equations  (3)  in  root  locus  form,  considering  My  as  a  variable 

quantity.  The  Laplace  transform  of  equations  (3)  is  taken,  and  the  charac¬ 
teristic  equation  is  calculated;  then,  the  characteristic  equation  is 
arranged  in  root  locus  form  as 


«U  (e  -  Vo.) 

- £ - l  (it) 

(»  -  -  [M-  +  m*u0  ]  s  -  n,u  ) 


Now,  the  characteristic  roots  of  the  quadratic  factor  in  the  denominator 
of  equation  (4)  have  been  determined  from  the  single-degree-of-freedom 


8 


tests  [equation  (l)]  and  are  shown  in  Figure  11  (^).  The  derivative  Xy 
is  then  specified  by  the  condition  that  the  root  locus  for  variable  My 

must  pass  through  the  experimentally  measured  data  (Q),  as  shown  in 
Figure  11. 

The  stability  derivative  My  is  determined  from  the  gain  [My  (g  -  X^Uq  )] 

required  for  the  calculated  roots  to  agree  with  the  experimentally  measured 
roots.  The  derivative  values  that  were  found  are  listed  in  Table  IV. 


Since  the  two-degree-of-freedom  experiments  were  conducted  with  various 
levels  of  rate  feedback,  a  verification  of  the  calculated  value  of  Xy  is 

possible.  It  is  assumed  that  fore  and  aft  differential  propeller  blade 
angle  produces  only  a  pitching  moment,  and  that  there  is  no  lag  between 
the  rate  signal  (0)  and  the  control  actuation  (A0p|TCH  )>  so  that  the  effect 


of  rate  feedback  may  be  included  as  an  increment  in  pitch  damping, 


The  increment  produced  by  the  rate  feedback  is  also  calculated.  The 

characteristic  equation  derived  from  equations  (3)  is  rearranged  in  a  form 

expressing  aMa  as  variable. 

0 


■  AMe  (s  -  Xy)  S 

[s  -  XyHs2  -  (M-  +  M*U0f)  s  -  MwU0f] 


-  =  -  l 

+  ^  (B  -  V0f) 


(5) 


The  roots  of  the  polynomial  in  the  denominator  of  this  expression,  i.e., 
the  poles  for  the  root  locus,  are  the  characteristic  dynamics  with  no 
feedback.  These  are  known  from  the  two-degree-of-freedom  analysis  with 
no  feedback.  The  zero  for  the  root  locus  based  on  equation  (5)  is  located 
on  the  real  axis  at  Xy.  Therefore,  the  value  of  Xy  is  verified  from  the 

condition  that  this  locus,  drawn  for  AM^  varying,  corresponding  to  various 

feedback  gains,  intersects  the  experimentally  determined  dynamics  for 
various  feedback  gains. 

The  increments  in  pitch  damping,  as  a  function  of  rate  feedback  gain, 
determined  by  this  procedure,  are  given  in  Table  IV.  The  value  of  Xy 

determined  by  this  approach  agreed  closely  at  all  duct  incidences  with  the 
value  determined  from  the  unstabilized  model  responses. 

Two-degree-of-freedom  motions,  consisting  of  pitching  and  vertical  ve¬ 
locity,  were  also  measured.  At  all  incidences  except  50°,  the  character 
of  this  motion  was  divergent;  thus,  no  analysis  was  attempted  because  of 
the  difficulty,  previously  mentioned,  of  making  quantitative  measurements 
of  divergent  motions.  At  50°  incidence,  the  characteristics  of  this  two- 
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degree-of-freedo*  motion  My  b«  used  to  determine  7V,  u  seen  fro* 
equation  (17)  In  Appendix  I,  with  xC(  and  «c<[  equal  to  aero,  u  rollout i 

4r  -  Vr  •  K\  *  *  0 

•  -  Vf  *  *  •  <v  tv 4  •  Vof  •  •  °  <«> 


Comparison  of  aquations  (6)  with  the  slngle-degree-of-freedo*  case 
[equation  (1)]  shout  that  there  la  one  additional  stability  derivative 
present  In  equations  (6);  that  derivative  is  Zy.  The  dovnvash  lag  derive* 

tlve  (M*)  nou  appears  separately  fra*  the  pitch  damping  derivative  (M*h 

It  is  therefore  possible  to  obtain  an  indication  of  the  size  of  this  de- 
rlvative  fro*  a  root  locus  dram  uith  Zy  as  a  variable  parsswter.  The 

root  locus  equation  developed  fro*  equations  (6)  takes  the  following  for*: 

•  2L  (•  -  »U)  s 

. — ..JL- - t. — - ,  .  .  i  (7 

■  (■■  -  fHj  *  H»U0f).  •  Vof> 

The  locations  of  the  two  poles,  determined  fro*  the  quadratic  factor  in 
the  denominator  of  equation  (7),  are  the  slngle-degree-of-freedo*  roots, 
equation  (1),  and  therefore  are  known.  The  pitch  dsaplng  is  determined 
from  the  condition  that  the  locus  of  roots  with  variable  Zy  Intersects 

the  experimental  points.  The  value  of  Zy  Is  calculated  fro*  the  gain 

required  for  coincidence  of  the  calculated  and  measured  roots.  At  the  one 
duct  incidence,  50°,  where  this  analysis  was  made,  indications  were  that 
was  negligible  compared  to  Mg. 

THREE  DEGREES  OF  FREEDOM 


The  three-degree-of- freedom  motions  were  analysed,  and  the  derivatives  Zy 
and  2^  were  calculated.  In  most  cases,  in  the  low-speed  flight  regime, 

these  derivatives  are  quite  small,  producing  only  saall  changes  between 
the  two-degree -of- freedom  (e-Uf)  motion  and  the  three-degree-of-freedo* 

(6-Uf-wf)  motion  (Figure  10). 

If  Xw  and  are  assumed  to  be  negligible,  then  the  three-degree-of- 

freedom  characteristic  equation  developed  from  equations  (15^,  in 
Appendix  I,  may  be  expressed  in  the  following  form: 
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4  '«■#*■*  /*W’W  • 


4e,u  ‘  ^  [(‘e.u'o  ’  'V  ^]  ■  0 


where  „  is  the  characteristic  equation  of  the  tvo-degree-of-freedcm 
aotion  and  „)  la  the  characteristic  equation  of  the  tvo-degree-of- 

vf  u  0 

freedom  motion  with  ■  0.  Fran  equation  (8)  it  may  be  noted  that  if  My 
is  equal  to  zero,  or  if  Zy  and  Zy  are  zero,  then  there  is  no  difference 

between  characteristic  dynamics  of  the  two-degree-of-freedom  motion 
involving  0  and  Uf  and  the  three -degree -of- freedom  motion  (e-Uf-wj. )  aside 

from  the  extra  root.  It  is  assumed  that  Xy  is  negligible,  and  this  is 

used  far  the  duct  incidences  in  which  there  are  measurable  differences 
between  the  two-  and  three-degree-of- freedom  motions.  The  polynomials 
A.  and  (d.  )  have  known  coefficients  from  the  two-degree-of-freedom 

9,u  e#u  o 

analysis.  Equation  (8)  Is  rearranged  in  root  locus  fora  as 
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Now  the  value  of  the  quantity  Zyg/Zy  will  determine  the  location  of  the 
zeros  of  equation  (9).  It  is  found  that  the  locus  of  roots  for  Zy  varying, 
at  a  constant  value  of  the  ratio  Zy/Zy ,  must  intersect  the  experimentally 
measured  values  of  the  frequency  and  damping.  The  value  of  Zy  is  calcu¬ 
lated  from  the  gain  required  for  coincidence  of  the  calculated  and  experi¬ 
mental  points.  Then,  Zu  is  determined  from  the  known  Zyg/Zy.  The  sta¬ 
bility  derivatives  determined  in  this  fashion  are  listed  in  Table  IV. 

This  procedure  was  generally  followed  at  all  duct  incidences,  to  evaluate 
the  stability  derivatives  of  the  vehicle,  with  minor  variations  as  noted 
below. 


Only  the  single-  and  two-degree-of-freedom  (e-uf)  motions  were  analysed. 

Experimental  results  showed,  as  would  be  expected  from  symmetry  consider¬ 
ations,  that  the  wf  motion  was  not  coupled  to  the  (e-uf)  aotion.  It  is 

not  possible  to  determine  Zy  in  hovering  because  of  the  nature  of  the 
tests  conducted. 

In  hover,  three  test  conditions  representing  different  blade  angle  and 
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propeller  rpm  settings  were  investigated.  Two  of  the  test  conditions 
utilize  different  combinations  of  blade  angle  and  rpm  to  produce  the  same 
total  thrust  (vertical  force  equal  to  the  weight  of  the  model):  8<7Bt 

equal  to  25.8  degrees,  rpm  equal  to  7000;  and  0t7B„  equal  to  29.2  degrees, 
rpm  equal  to  6400.  The  third  test  condition  uses  another  combination  of 
blade  angle  and  rpm,  resulting  in  a  lower  total  thrust  than  previous  cases: 
8.76*  equal  to  25.8  degrees,  rpm  equal  to  6400.  This  combination  of  blade 
angle  and  rpm  produces  a  hover  thrust  of  43.1  pounds,  corresponding  to  a 
scaled  gross  weight  of  14,000  pounds  for  the  Bell  X-22A.  Note  that  the 
weight  of  the  model  is  51 .5  pounds  in  all  cases,  so  that  only  two-degree- 
of- freedom  (0-uf)  motions  were  examined  in  the  low  thrust  case. 

Dimensional  analysis  can  be  used  to  demonstrate  that  for  the  two  test 
conditions  at  the  same  propeller  blade  angle,  if  the  time  scale  of  the 
dynamics  is  nondimens ionali zed  by  the  rpm,  then  the  nondimens ional  frequency 
and  damping  should  be  Independent  of  the  rpm.  This  comparison  of  the  two 
test  cases  is  shown  in  Figure  10.  The  spread  in  the  points  is  considered 
to  be  within  the  accuracy  with  Which  frequency  and  damping  can  be  evaluated 
from  the  highly  unstable  dynamics  of  this  hovering  motion. 

If  the  dynamic  stability  of  the  vehicle  depends  only  on  the  geometric 
configuration  of  the  vehicle  and  the  duct  exit  velocity,  then  there  should 
be  no  difference  between  the  dynamic  characteristics  at  the  same  total 
thrust  level  produced  by  different  combinations  of  blade  angle  and  rpm. 

There  are  measurable  differences,  although  not  large,  between  these  two 
test  conditions.  The  results  indicate  that  it  is  desirable  to  conduct 
dynamic  stability  tests  at  the  proper  blade  angle. 

V600 

The  proced  ure  described  above  for  id  -  90°  was  followed.  No  measurable 
difference  existed  between  the  two-degree- of- freedom  motion  (0-Uf)  and  the 
three -degree -of -freedom  motion  (0-Uf-vf );  thus,  Zy  and  were  not 

determined.  That  is,  similar  to  the  hovering  case,  there  is  only  weak 
coupling  between  the  (0-Uf)  degrees  of  freedom  and  the  w^  degree  of  freedom. 

It  is  expected  that  2y  has  a  value  in  this  flight  condition  as  it  has  in 
hover,  but  it  cannot  be  determined  from  the  experiments  conducted. 

*d  ■  70° 

The  analysis  at  this  trim  condition  is  discussed  in  detail  in  Appendix  II, 
and  the  transient  response  data  are  shown  in  Figures  17  through  28.  The 
transient  response  of  the  modal  in  two  degress  of  freedom  (0-U^),  shown  in 

Figure  18,  exhibited  a  rate  of  growth  of  the  dominant  unstable  oscillation 
(with  no  rate  feedback)  which  was  so  rapid  that  it  was  difficult  to  measure 
the  amplitude  ratio  from  the  time  history,  which  extends,  at  most,  for  one 
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cycle.  Thus,  the  experimental  results  for  one  particular  rate  feedback 
setting  (K»  ■  0.041+  sec )  In  two  degrees  of  freedom  were  analyzed.  It  was 

assumed  that  the  Increment  in  pitch  damping  produced  by  this  feedback  gain 
was  equal  to  that  produced  by  the  same  gain  setting  at  80°  duct  incidence. 
This  known  Increment  in  pitch  damping  is  added  to  the  value  determined 
from  the  single -degree -of- freedom  tests.  Then  the  analysis  proceeds  as 
described.  The  dynamics  of  the  vehicle  with  no  feedback  were  then  calcu¬ 
lated  from  the  resulting  derivatives.  The  resulting  characteristic  roots 
agreed  with  the  information  that  could  be  determined  from  measured 
transients.  The  increments  in  damping  corresponding  to  other  feedback 
settings  agreed  closely  with  those  in  the  80°  case,  as  shown  in  Table  IV. 
This  confirms  the  assumption,  inherent  in  this  approach,  that  the  moment 
produced  per  degree  of  differential  fore  and  aft  blade  angle  at  id  =  80° 
is  equal  to  that  at  ld  •  70° . 

Id  ~  60° 

It  is  necessary  to  add  certain  known  terms  to  the  equations  of  motion  at 
id  ■  60°  to  account  for  the  displacement  of  the  eg  of  the  model,  with 

respect  to  the  pivot  axis,  as  explained  in  Appendix  I.  Values  of  these 
terms,  Which  were  added  to  the  equations  of  motion  before  proceeding  with 
the  analysis,  are  listed  in  Table  IV.  No  alteration  in  the  analysis 
procedure  is  required,  since  only  known  terms  are  added  to  the  equations 
of  motion. 

In  addition,  a  complication  associated  with  the  model  control  system  was 
encountered  at  this  flight  condition.  When  the  carriage  was  commanded  by 
the  model  to  accelerate  to  follow  the  model  motion  at  these  relatively  high 
trim  speeds,  the  increased  current  drawn  by  the  carriage  drive  motor  caused 
a  noticeable  drop  in  line  voltage.  This  line-voltage  drop  affected  a  power 
supply  in  the  model  control  system,  producing  a  propeller  blade  angle 
change  of  equal  value  on  all  four  oropellers  approximately  proportional  to 
carriage  (model)  acceleration.  The  configuration  of  the  vehicle  is  such 
that  total  propeller  blade  angle  change  causes  pitching  moments  (Reference 
3),  thus  providing  an  apparent  pitching  moment  variation  with  horizontal 
acceleration.  It  was  assumed  therefore  that  this  effect  could  be  accounted 
for  by  adding  the  unknown  derivative  to  the  analysis.  It  is  then 

possible  to  analyze  the  data,  remove  this  effect,  and  determine  the  sta¬ 
bility  derivatives  by  an  analysis  similar  to  that  described  above. 

Modifications  to  the  model  control  system  were  made  such  that  this  coupling 
phenomenon  was  eliminated  in  the  50°  duct  incidence  tests  and  all  future 
tests. 

*d  -  5°° 

The  50°  case  also  required  additional  terms  in  the  equations  of  motion  to 
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include  the  effect  of  the  displacement  model  eg  from  the  pivot  axis.  Also 
as  mentioned  earlier,  no  mechanical  springs  were  used  in  the  single -degree 
of-freedom  tests  since  the  angle -of -attack  stability  derivative  was  nega¬ 
tive  (MwUo  <  0).  Otherwise,  the  analysis  proceeds  as  described. 

Table  IV  lists  the  stability  derivatives  found  in  model  scale  based  on 
the  inertia  characteristics  listed  in  Table  III.  The  model  scale  deriva¬ 
tives  are  also  shown  graphically  in  Figure  12. 
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STABILITY  DERIVATIVES  OF  THE  FULL-SCALE  AIRCRAFT 


The  stability  derivatives  determined  for  the  model  can  be  interpreted  in 
terms  of  full-scale  aircraft  characteristics 

It  is,  of  course,  necessary  to  assume  that  there  are  no  scale  effects  to 
make  this  interpretation.  Comparison  of  the  lift,  drag,  and  pitching 
moments  on  an  isolated  duct  of  the  model  and  the  full-scale  aircraft  re¬ 
vealed  that  there  were  no  appreciable  scale  effects  on  the  model  duct  in 
the  flight  conditions  of  interest  (Reference  1).  Therefore,  important 
scale  effects  are  not  expected  to  be  present  in  the  results. 

Rather  than  present  nondimensional  derivatives,  it  is  considered  to  be 
more  convenient  and  conventional  to  present  dimensional  force  derivatives 
divided  by  the  mass  of  the  aircraft,  and  pitching  moment  derivatives 
divided  by  the  inertia.  The  full-scale  derivatives  are  based  on  the 
moment  of  inertia  of  the  Bell  X-22A,  as  given  in  Reference  3« 

The  dimensional  full-scale  derivatives  are  listed  in  Table  VI  and  are 
shown  graphically  in  Figure  13  for  the  full-scale  aircraft,  using  the 
altitude  gross  weight  equivalence  discussed  in  Appendix  III.  The  relation¬ 
ship  that  applies  for  the  test  program  is  shown  in  Figure  l4,  and  the 
conversion  factors  that  result  are  given  in  Table  V.  Thus,  for  example, 
the  derivatives  presented  correspond  to  the  X-22A  flying  at  a  gross  weight 
of  16,800  pounds  at  sea  level  or  a  gross  weight  of  li+,000  pounds  at  a 
density  altitude  of  6000  feet. 

The  stability  derivatives  show  the  following  trends. 

THE  SPEED  STABILITY  (My) 


This  derivative  is  large  and  positive  at  the  three  highest  duct  incidences 
tested;  it  decreases  considerably  at  duct  incidences  of  60°  and  50°.  The 
large  value  near  hovering  and  in  low- speed  flight  is  the  primary  source  of 
the  oscillatory  instability  present  in  the  data;  also  it  indicates  that 
the  vehicle  will  exhibit  an  appreciable  sensitivity  tD  horizontal  gusts. 

THE  ANGLE -OF -ATTACK  STABILITY  (MwUq  ) 


This  derivative  is  fairly  large  and  positive  (unstable)  at  all  but  the 
lowest  duct  incidence  tested.  At  a  duct  incidence  of  50°,  the  derivative 
is  negative.  The  unstable  value  of  the  derivative  contributes  to  the 
instabilities  of  the  motion  at  the  higher  duct  incidences,  and  the  change 
in  sign  is  the  primary  contributor  causing  the  stable  motion  at  a  duct 
incidence  of  50°.  The  trend  in  this  derivative  is  typical  of  V/STOL 
aircraft  at  low  speeds. 
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THE  PITCH  DAMPING  (M^) 

This  derivative  is  comparatively  small  and  negative  and  generally  increases 
with  decreasing  duct  Incidence.  It  is  not  clear  at  this  time  why  the  trend 
in  this  derivative  does  not  follow  a  smooth  curve.  The  low  values  Indicate 
that  the  full-scale  vehicle  would  require  damping  augmentation  to  provide 
satisfactory  handling  qualities  at  low  speeds.  The  small  value  of  this 
derivative  may  be  seen  by  noting  the  large  increments  in  damping  (roughly 
a  factor  of  10  at  duct  incidences  of  80°  and  70°)  required  to  make  the 
transient  motion  neutrally  stable. 

THE  DOWNWASH  LAO  (M$) 


All  indications  from  the  data  are  that  this  derivative  is  small  compared 
to  the  pitch  damping. 

THE  RATE  OF  CHANGE  OF  HORIZONTAL  FORCE  WITH  HORIZONTAL  VELOCITY  (\) 


This  derivative  is  large  and  negative  at  low  speeds.  The  primary  source 
of  this  der  vative  is  the  momentum  drag  of  the  ducts. 

THE  LIFT  CURVE  SLOPE  (Zw) 

This  derivative  is  small  and  increases  with  speed.  The  value  of  the 
derivative  was  not  determined  in  hovering.  The  values  at  low  speeds  must 
be  considered  as  approximate,  since  there  is  only  weak  coupling  between 
the  two-degree-of -freedom  (g-Uf)  motion  and  the  three -degree -of- freedom 
(8-Uf-Vf)  motion.  This  makes  it  difficult  to  determine  Z^  accurately. 

THE  RATE  OF  CHANGE  OF  VERTICAL  FORCE  WITH  HORIZONTAL  VELOCITY  (Z^) 


This  derivative  is  small  and  of  normal  sign  (negative)  in  the  cases  where 
it  was  evaluated. 

THE  RATE  OF  CHANGE  0*  HORIZONTAL  FORCE  WITH  VERTICAL  VELOCITY  (X^.) 


This  derivative  is  normally  small  and  was  not  determined  from  the  experi¬ 
ments.  It  was  assumed  to  be  negligible  in  the  analysis. 
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CONCLUSIONS  AND  RECOMMENDATION 


CONCLUSIONS 

1.  The  pitching  moment  derivatives  of  this  quad  configurat  on  v/STOL 
model  exhibit  the  following  characteristics: 

a.  The  speed  stability  (My)  is  large  and  positive  (statically  stable) 
at  high  duct  incidences. 

b.  The  angle-of-attack  stability  (M^Uq  )  is  positive  (unstable)  at 
all  but  the  lowest  duct  incidence  tested  (50°). 

c.  The  damping  in  pitch  (M*)  is  stable  but  small  in  hover  and 
increases  with  decreasing  duct  incidence. 

2.  For  analysis  of  the  data,  the  inclusion  of  test  conditions  with  pitch 
rate  feedback  was  valuable,  particularly  in  the  experiments  where  the 
basic  model  was  highly  unstable. 

RECOMMENDATION 


It  is  recommended  that  an  effort  be  made  to  correlate  the  results  obtained 
herein  in  the  form  of  stability  derivatives  with  full-scale  flight-test 
data  on  similar  configurations. 


TABLE  I.  SUMMARY  OF  LONGITUDINAL  TEST  CONDITIONS  AND  ME 


Continued 


Duct  Average  Differential  Propeller  Trim  Degrees  Stability  Roots  of 

Incidence  Propeller  Collective  Speed  Velocity  of  Augmentation  Oscillatory 

Pitch  Pitch  Freedom  Motions 
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TABLE  I  -  Concluded 


TAHLE  II.  SCALE  FACTORS  KOh  !'V!.AM1C  MOUKL  SIMILARITY 


Multiply  full-scale  property  by  scale  factor 

to  obtain  modsl  property. 

For  Xt  ■  0.1,i‘>1 

Linear  dimension 

\ 

0.1*i*>3 

Area 

X  9 

At, 

2.112  x  10-* 

Volume,  mass,  force 

K* 

3.071  x  10  •» 

Moment 

\4 

x  if* 

Moment  of  inertia 

<»/87  x  10 -* 

Linear  velocity 

»  O.S 

Al 

0. 3812 

Linear  acceleration 

\° 

1.000 

Angular  velocity 

%  ■  0*5 

*.&?3 

Angular  acceleration 

K 

G.1V>3 

Time 

1  02 

A», 

\  38L’ 

Frequency 

\  -0.S 

Al 

2/23 

Reynolds  number 

v  1.6 

'l 

•  .''-l  x  1C’5 

Mach  number 

\  0.S 

Al 

0.3812 

where  X,  ■ 

me  del  linear  dimension 

full-scale  linear 

dimension 

TA3LF  III.  WOOCL  OEOMmiC  AJC  ZSSariA  CH 


50  1  ^.70 1  20.10 1  2.^6  |  0.67 

*81ml«~^*cr— -of-fp— do«  runs  only. 

**NultipI«  d«fr«*-of- fpMdoa  runs  only. 


TABU  IV.  MOONL  STANII.IT*  UNIVATIVKS  FOK  TKIiKED  LKVKU  KLIQHT 


Duct  Incidence 
(deg) 

90 

80 

70 

60 

50 

f  Model  Trio  Velocities  end  Aerodynamic  Derlvetlves 

0 

10 

22 

28 

36 

-.53 

-.**5 

-.70 

-  .60 

-.50 

Assessed  Negligible  in  all  Cases 

Neg 

Neg 

*.19 

-.655 

-.830 

2E5M 

NA 

NA 

-.11* 

-.675 

-.895 

3IE2823I 

♦.307 

♦.1*21 

♦.1*81* 

+.157 

+.154 

NA 

♦.255 

+.297 

+.110 

-.160 

-.1*1* 

•  .1*8 

-.955 

-.51* 

-.864 

Additional  Stability  Derivatives  Due  to  eg  Offset 

mm 

0 

0 

0 

-.0304 

-.034 

an 

1  0 

0 

0 

+.125 

+0.13 

mm 

0 

c 

0 

-.98 

-1.09 

Additional  Stability  Derivatives  Due  to  Stsbility  Augmentation 

&Mq»  1/*«c 

NT 

NT 

NT 

-2.59 

NT 

g 

NT 

-*♦.37 

-1*.  35 

NT 

NT 

NT 

im 

-8.15 

NT 

NT 

BB 

NT 

-11.9 

-11.02 

NT 

NT 

an 

0 

0 

0 

+.075 

0 

Abbreviations:  NA  Not  Available 

Neg  Negligible 

NT  Not  Tested 

Test  Conditions  given  in  Table  I 

Inertia  Characteristics  given  in.Table  III 

Trie  Conditions  alven  InTimre  8  ------ 

TABLE  V.  INTERPRETATION  OF  FORCES,  MOMENTS,  AND  VELOCITIES  AT 
OTHER  GROSS  WEIGHTS 


Altitude 

Gross  Weight 

Velocity 

Gross  Weight 

Forces 

V 

V 

Moments 

V 

Aw 

Velocities,  angular 
and  linear 

1 

Aw*5 

Air  density 

Aw 

1 

Angles 

1 

1 

wo 

desired  gross  wei$ 

;ht 

'V  Wc  gross  weight  determined  by  dynamic  scaling 

To  determine  aerodynamic  quantities  at  other  gross  weights,  multiply 
dynamic  scaling  results  by  the  above  quantities. 

NOTE:  Use  of  the  first  column  results  in  no  change  in  dynamic  stability 
characteristics.  Use  of  the  second  column  results  in  changes  in 
dynamics . 
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TABLE  VI.  FULL-SCALE  STABILITY  DERIVATIVES  FOR  TRIMMED  LEVEL  FLIGHT 


Duct  Incidence 
(deg) 


U0  ,  ft/sec 

0 

26.2 

57.5 

73.4 

94.2 

l/ sec 

-.20 

-.17 

-.27 

-.23 

-.19 

V  1/sec 

Assumed  Negligible 

Zu>  1/sec 

Neg 

Neg 

-.073 

-.25 

-.32 

l/sec 

NA 

NA 

-.054 

-.26 

-.34 

Mu,  l/ft-sec 

+.017 

+.026 

+  .028 

+.011 

+.011 

M^,  l/ft-sec 

NA 

+  .015 

+.017 

+.0077 

-.0112 

M»,  l/sec 

-.17 

-.20 

-.40 

VO 

CO 

• 

1 

-.42 

Abbreviations:  NA  Not  Available 

Neg  Negligible 

See  Figure  14  for  density  altitude/ gross  weight  correspondence 
Radius  of  gyration  ky  =  8.5  ft. 

'  • 
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Figure  1.  .‘rxnceton  Dynamic  Model  Track  -  Model  Mounted 
on  Longitudinal  nynamic  Testing  Apparatus. 
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Figure  2.  0.145  Scale  Quad  Configuration  Dynamic  Model. 


Figure  4.  Location  of  Model  Reference  Stations  and  eg 


RADIAL  STATION,  r/R 


Figure  5.  Geometric  Characteristics  of  Ihree-Blmded 
Model  Propellers. 
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CHORD-  DIAMETER 
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DUCT  INCIDENCE,  ^,d«g  DUCT  INCIDENCE,  ld,d«g 


Figure  13.  Full-Scale  Stability  Derivatives  and  Trim  Condition  for 
Altitude/Gross  Weight  Correspondence  Shown  in  Figure  1^ 
(ky  =  8.5  ft). 
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APPENDIX  I 
FIXATIONS  01  WflON 


Linearized  equations  of  motion,  applicable  to  the  analysis  of  various  ex¬ 
perimentally  measured  responses,  are  presented  In  this  appendix. 

The  longitudinal  equations  of  motion  that  describe  the  small  perturbation 
motion  of  an  aircraft  frn*r.  initially  level  flight,  using  a  stability  axis 
system  (nofirenji  C),  arc. 

U  -  XyU  -  XyW  +  g  0  «*  0 

w  -  ZyW  -  ZyU  -  U0  0  =  0 

MyW  +  M^w  +  N^u  +  Mg0  -0=0  (10) 

IVo  derivatives  X*  and  Z*  that  are  usually  small  are  neglected. 

0  0 

Since  all  of  the  transient  responses  were  measured  and  are  presented  in 
terms  of  space-fixed  variables,  it  is  convenient  to  transform  equations 
(10)  to  a  space-fixed  system  (Figure  15),  with  the  Xf  axis  parallel  to  the 
horizon,  by  -/he  following  transformations: 

u  =  uf  -  w0f>  e 

w  =  wf  +  Uc^  0  (11) 

where  Wn  is  equal  to  zero  from  the  condition  of  initially  level  flight. 

'f 

Substituting  the  relationships  of  equation  (11)  into  equations  (10),  the 
following  equations  result: 

uf  -  Xuuf  -  XyWf  +  (g  -  X^)  0  =  0 

*f  “  Vf  -  Vf  •  V^f  6  -  0 

e  -  (*$  ♦  6  -  tyj0f  e  -  ^  -  M**f  -  Mjtf  =  o  (12) 


•*  v*,  a  >,}.'?*■  ■  ■ 

•  rf '  >  1  *  ’  '  . 


» f  -tw/if  **  iii  n 


.  tfCJKJK 


iMMtft  naitw  tmm  wr ■irtf 


!*.-.  i  ufcfe  n  certain  features  of  the  model  and  the  apparatus,  three  modifi¬ 
cations  to  these  equation.^  are  necessary  such  that  they  will  apply  tc  all 
test,  conditions. 

1.  There  are  two  linkages  required  to  attacn  tne  model  Lw  tiiw  wOrv* 
transducers  and  mounting  system  used  for  this  type  of  testing. 
These  supports  provide  the  horizontal  and  vertical  translational 
degrees  of  freedom  and  contribute  additional  masses  (m^  and  niy) 

that  "fly"  along  with  the  model  and,  therefore,  must  be  acceler¬ 
ated  by  the  model.  The  two  linkages  art  relatively  light  in 
ve^ht  composed  to  the  "flying"  weight  s  f  the  model  but  never- 
i-hrti*ipg  c-ho’ild  be  amounted  for  Dy  additional  mass  terms  in  the 
equations  of  motion.  Generally,  the  arrangement  and  weights  of 
these  two  supports  are  such  that  the  mass  accelerated  by  the 
model  iii  the  horizontal  direction  is  larger  than  that  accelerated 
in  the  vertical  direction.  If  mp  is  the  total  mass  of  the  model 

resting  on  the  pivot  axis  (Figure  16),  then  the  total  lifted  mass 
of  the  model  m  when  "flying"  is  equal  to  mp  plus  the  mass  of  the 

vertical  link  or  m  =  nip  +  itiy.  Similarly,  the  total  accelerated 
macs  in  the  horizontal  direction  n^.  is  equal  to  mp  +  niy.  +  or 
m  +  mh.  This  d„ uamic  model  mount  characteristic  requires  the 

..  modification  of  all  terms  in  the  horizontal  force  equat*. ~n,  except 
the  acceleration  term,  by  a  mass  ratio  defined  as  m/mj.  and  equal 

v  to  0.936  in  value. 

?..  certain  of  the  test  conditions,  as  indicated  in  Table  Ilx,  tug 
ct  :f  the  model  was  not  located  at  the  pivot  axis  of  the  model. 
Equations  (12)  may  be  considered  to  be  written  about  the  pitch 
piv.u'  axis  of  the  model,  which  represents  the  full-scale  eg 
position  about  which  the  derivatives  are  determined.  Additional 
terms  a’-e  necessary  in  the  equations  motion  to  account  for  the 
displacement  0/  the  model's  eg.  These  are: 


W  Z 
P  eg 


where  Mp  and  Wp  are  re*,  actively  the  pi  voting  mass  and  pivoting 
weight  cf  the  model. 


3.  In  certain  of  the  tests  (tingle  degree  of  freedom  only),  a 

mechanical  apring  vat  added  about  tho  model  pitch  axia  to  provide 
a  restoring  moment  which  produces  an  oacillatory  motion  of  tho 
model.  In  theae  experiments  the  following  term  should  be  added: 


(1M 


In  the  experiments  Where  a  apring  was  employed,  the  value  of  the 
apring  constant,  k~  ,  is  as  given  in  Table  III. 
em 

Adding  the  necessary  terms  to  account  for  these  three  effects,  the  com¬ 
plete  equations  of  motion  that  apply  to  the  measured  transients  obtained 
in  this  facility  are: 


-  s;  V»f  -  a;  Vf  ♦  ^  <*  -  ¥»f)  8  •  ° 

Wf  -  Zy.Wj  -  ZyUj.  -  0  *=  0 


This  set  of  equations  would  apply  for  the  three-degree-of-freedom  tests  if 
the  k  terms  were  removed. 

e_ 

m 

For  the  restricted  degree  of  freedom  tests,  the  following  reduced  sets  of 
equations  apply. 

1.  In  two  degrees  of  freedom,  with  kfl  =  0: 

0,  uf  (wf  -  0) 

*r  '  ^  Vr  *  ^  <«  -  V.f)  8  ■  0 
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*  -  <1(1  *  Vof>  4  ■  Vof  •  ♦  ^  <if  -  M„«f  .  0 

(16) 


e,  vf  (ut.  -  o) 

wf  "  ^wwf  "  8  “  0 


V  -  *  MjU0f)  e  -  tyj0f  e 


+ 


W-^«e 

xy 


o 


(17) 


2.  In  the  single-degree-of- freedom  experiments,  with  the  mechanical 
spring  and  uf  *  0,  wf  -  0,  the  equation  that  applies  Is 

e  -  (^  +  M*u0f)  i  +(i»  -  iy>0f  +  e  -  o  (16) 


3.  In  the  experiments  where  feedback  is  used,  a  term  MA(I  dp 

flPMTCM  f,TCH 

should  be  added  to  the  right-hand  side  of  the  pitching  moment 
equation,  and  then  the  equation  governing  A0  is 


■  k4  4  (19) 


By  substitution  of  these  expressions  into  the  pitching  moment 
equation,  an  effective  pitch  damping  is  obtained 


pi  TCH 
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(20) 


SPACE-FIXED  AXIS 


U,  =H,f+u, 


STABILITY  AXIS 


(body  fixed;  initially  aligned  with  tr«»«tream 
velocity  at  forward  epeatfe  or  with  hjrizon  in  hover) 


Figure  15.  Definitions  of  Space-Fixed  and  Stability  Axis  Systems 

'Variables  Are  Shown  in  Their  Positive  S°nse). 
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MODEL  AND  LINKAGE  MASS  ARRANGEMENT 


LifIM  mm  ■  >  Mp*av 

T»M  NniaMaiii'  a,  «af  *  aT«  a^|«a»ak 


MODEL  Cl  -  PIVOT  AXIS  REFERENCE  SYSTEM 


l 

Figure  16.  Model  and  Link  Mass  Arrangement  and  Reference  System 
for  Model  eg  and  Pivot  Axis. 
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APPENPIX  II 
NUMERICAL  EXAMPLE 


In  order  to  illustrate  in  detail  the  method  used  to  analyze  the  data,  the 
calculations  at  a  duct  incidence  of  70°  are  presented  in  this  appendix. 
Data  at  this  duct  incidence  from  Reference  2  are  repeated  here  in  Figures 
17  through  28. 

First,  the  s ingle -degree- of- freedom  data  were  analyzed.  From  Appendix  I, 
the  equation  of  motion  is 

6  "  (M§  +  M*U0f)  «  -  *fUof  ♦  6  -  0  (21) 


The 


mechanical  spring  constant, 


is  determined  by  calibration  to  be 


18.1  foot-pounds  per  radian;  Wp  zcg  is  found  to  be  I.89  foot-pounds  per 

radian  and  arises  from  the  fact  that  the  eg  of  the  model  is  below  the 
pitching  axis  for  this  single-degree-of-fi eedom  test.  With  the  model's 
motor  off  and  the  forward  velocity  equal  to  zero,  measurements  were  made 
to  determine  the  moment  of  inertia.  The  aerodynamic  damping  is  assumed 
to  be  negligible,  and  the  equation  of  motion  in  this  case  is 


e  ♦ 


iv  v  iv  ) 


e 


(22) 


where  is  mechanical  lamping  due  to  bearing  friction.  The  measured 
roots,  determined  from  Figure  17,  are 


a  ■  -  0.029  *  2.661  per  second 


Thus,  the  characteristic  equation  is 


s*  ♦  0.0>8s  ♦  7.15  -  0 


(23) 


Then  the  moment  of  inertia,  Iy,  is  fourd  r'r: 


vfi  ♦ 

oB  p  eg 


7.15  per  second  squared 


*  »*fcr"*?*^«**r*‘' 


ft r-  r^'j'J^lV^*™******^^ 


and  is  therefore  equal  to  2.80  slug-feet  squared,  and  the  mechanical 
damping,  Cm/ly,  is  equal  to  0.058  per  second.  Balance  weights  present  in 
the  model -motor-off  tests  produced  a  moment  of  inertia  increment  of  0.65 
slug-feet  squared.  The  moment  of  inertia  of  the  model  for  the  tests  at 
forward  speed,  with  propellers  running,  is  therefore  2.15  slug-feet  squared. 

At  the  trim  speed  for  duct  incidence  equal  to  70°  (22  feet  per  seconi)  and 
rpm  equal  to  6780,  the  single-degree-of-freedom  roots  measured  from  figure 
17  are 


s  =  -  0.57  ±  1.57i  per  second 


The  characteristic  equation  is 

s3  +  l.l4s  +  2.785  =  0  (24) 

Thus,  by  comparison  with  equation  (18), 

W  2  ^0 

c8  +  -  M^Uq  *  2.785  per  second  squared 

*y  *y  f 

-  (^  +  M*  +  M*U0  )  =  l.l4  per  second 

Iy  o  f 

Substituting  for  Wp,  zcg,  k  ,  Iy,  and  Ca, 

“m 

=6.62  per  second  squared 
-  (M-  +  M^U0  )  =  1.07  per  second 


In  all  trim  conditions,  as  in  this  example  case,  the  mechanical  damping 
was  small  in  comparison  to  the  aerodynamic  damping. 

The  pendulous  and  mechanical  spring  terms  are  now  removed  from  the  equation 
of  motion  in  order  to  proceed  to  the  two-degree-of- freedom  case.  The  re¬ 
sulting  single-degree-of-freedom  equation  is 

s8  +  l.l4s  -  6.62  =  0  (25) 

The  moment  of  inertia  in  the  two-  and  three-degree-of-freedom  cases  was 
slightly  different  from  the  single-degree-of-freedom  case  due  to  balance 


-<VtT-W 
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weights  and  is  equal  to  2.23  slug-feet  squared.  With  this  modification, 
the  above  equation  becomes 

s3  +  1.09s  -  6.39  =  0  (26) 

Now,  in  the  two-degree-of-freedom  case,  the  motion  of  the  model  was  highly 
unstable,  as  shown  in  Figure  18,  making  it  difficult  to  measure  the  period 
and  damping  when  there  was  no  feedback.  The  two-degree-of-freedom  case 
with  a  feedback  gain  K*  equal  to  O.ObU  second,  shown  in  Figure  22,  was 

analyzed,  assuming  that  the  damping  increment  produced  by  feedback  (see 
Appendix  I)  was  equal  to  that  at  the  same  feedback  gain  at  a  duct  incidence 
of  80°.  The  increment  at  a  duct  incidence  of  80°  had  been  determined 
previously.  From  Appendix  I,  this  damping  increment  AMg  is  equal  to 

K*M  and  therefore  is  equivalent  to  assuming  that  the  control 

"  I  TC  M 

effectiveness  is  the  same  at  80°  as  at  70°. 

AMg  was  determined  to  be  equal  to  -  8.15  per  second  at  80°  duct  incidence. 

This  value  of  damping  is  added  directly  to  equation  (26)  to  yield 

s2  +  9.21+s  -  6.39  «  0  (27) 


The  roots  of  this  equation  are 

Sj  =  -  9*89  per  second 
s2=  +  C.65  per  second 

Now  we  proceed  to  analyze  the  two-degree-of-freedom  motion.  Note  that  at 
this  duct  incidence  the  model  eg  was  coincident  with  the  pivot  axis. 

The  equations  of  motion  are,  from  Appendix  I,  (with  x„_  and  z  equal  to 
\  c8  c8 
zero ; , 

+  ^  -  Wf)  e  -  o 

-  Muuf  +  e  -  (m.  +  h*u )  e  -  iyj0  b  *  0  (28) 


The  notation  Mg 


implies  that  the  damping  includes  the  effect  of  feedback. 


m/m,  is  the  factor  due  to  the  model  mounting  linkage  and  is  equal  to  0.936. 

w 


Taking  the  Laplace  transform  of  the  above  equations,  and  placing  the 
characteristic  equation  in  root  locus  form,  considering  Mll  as  the  unknown 
parameter,  we  obtain 
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(29) 


0.936  ^  (g  -  X„U0f) 

(s  -  0.936  Xu)  (*  -  (Mg  +  M^u’f)  s  -  H,U0  ) 


The  quadratic  factor  in  the  denominator  is  that  determined  from  the 
single -degree- of -freedom  tests  [equation  (l6)]  and  has  roots  equal  to 
-9.89  per  second  and  +  0.65  per  second. 

Placing  these  known  quantities  in  the  above  equation,  we  obtain 

0.936  ^  (g  -  VJ0  ) 

- - -  3-1  [30) 

(s  -  0.936  Xu)(s  ♦  9.09)(s  -  0.65) 

The  measured  oscillatory  dynamics  for  this  two-degree-of-freedom  case  are 
(at  *  0.044  second) 

s  «  +  0.05  ±  l.OUi  per  second 

The  known  poles  (  ^ ^  from  the  single-degree-of-freeaom  analysis  and  the 
experimentally  measured  roots  ( O  )  are  shown  on  the  complex  plane  in 
Figure  29. 

Now,  there,  is  one  unknown  pole  located  at  0.936  ^  (©)•  The  location  on 

the  real  axis  is  determined  from  the  condition  frco  equation  (jO):  for 
positive  My,  the  sum  of  the  angle  contributions  from  the  Doles  (^,0  ) 

at  tlie  roots  (  O  )  must  be  equal  to  l80°.  From  Figure  29,  the  angle 
contributions  from  the  two  known  poles  are 


s  ■  ♦  0.65  per  second  Oj  *  117° 

s  ■  -  9.89  per  second  aa  *  6° 

The  unknown  pole  (0.936  Xy),  therefore,  must  produce  an  angle  contribution 

of  57°  to  make  these  three  angles  add  up  to  l80°.  It  is  located  at  -  0.6 
per  second,  as  shewn  on  Figure  29.  This  calculation  yields  a  value  of 
^  equal  to  -  0.6k  per  second.  Now  the  root  locus  for  varying  My  is 

drawn  as  shown  on  Figure  29,  and  the  gain  at  the  experimental  tvo-degree- 
of  freedom  roots  ( O  )  is  calculated.  This  will  determine  the  product 
My  (g  -  X^Uq  ).  The  gain  calculation  is 

,  %  (1.17)(1.23)(10) 

Hy  (g  -  XyUjj  )  «  - G.9y&  "  15.40  per  second  cubed 
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The  isolated  du<  t  data  of  Reference  1  indicate  that  Xyllg  at  this  flight 
condition  is  negligible  compared  to  g  and,  therefore, 

Mjj  ■  0.1*78  per  foot-second 


With  this  information  we  can  calculate  the  real  root  that  corresponds  to 
the  two-dagree-of- freedom  notion  from  the  characteristic  equation.  Taking 
the  calculated  values  of  My  and  Xy  and  placing  them  into  equations  (28), 
the  characteristic  equation  is  calculated 

s*  +  9.92  «*  -  0.52s  +  11.73  =  0  (31) 

The  roots  of  equation  (31)  are 

s,  ■  -  10.02  per  second 

sJf,  ■  +  0.05  ±  1.04i  per  second 

How  we  return  to  the  equations  of  motion  (28)  and  rearrange  the  character¬ 
istic  equation  in  root  locus  form,  considering  ^  as  a  variable  parameter 

to  determine  the  agreement  among  the  derivatives  found  at  a  K  •  =  0.44 

w 

second  and  the  other  two-degree -of- freedom  cases  with  different  levels  of 
rate  feedback.  In  addition,  the  data  can  be  extrapolated  to  calculate  the 
dynamics  of  the  vehicle  with  no  feedback.  The  measured  dynamics  at  other 
feedback  gains  from  Figures  21  and  23  are 

■  0.030  second  s  =  +  0.31  ±  1.251  per  second  (  fe)  ) 

Kj  ■  0.060  second  s  **  -  0.01  ±  0.93i  per  second  (^*) ) 


The  equation  for  the  root  locus  diagram  for  variable  M*  is 

U 


(•  dM^)  »(■  •  0.936  Xy) 

(s  -  0.936  V  (*  -  (^  ♦  M*u0f)  S  -  jyj0f))  +  ^  (g  -  vjGf) 


1 

(32) 


All  of  the  stability  derivatives  in  this  expression  have  been  determined 
except  AMg  and  may  be  substituted  giving 

(-  AMg)  s(s  +  0.60) 

(s  +  10.02)(s  -  0.05  +  1.04i)(s  -  0.05  -  1.04i) 


(33) 


The  root  locus  for  variable  AMg  may  now  be  sketched  as  shown  in  Figure  29. 

Note  that  this  locus  provides  a  verification  of  the  previously  calculated 
value  of  Xu,  since  it  must  pass  through  the  other  experimental  points  for 

different  feedback  gains.  The  0°  locus  ( - )  shows  the  trend  for  de¬ 

creasing  feedback  gain,  and  the  l80°  locus  ( — )  shows  the  trend  for  in¬ 
creasing  feedback  gain.  The  root  locus  passes^through  the  other  two 
experimentally  measured  characteristic  roots  (0  Q  )  and  verifies  the 
value  of  X^  Now,  the  increment  in  damping  provided  at  the  experimental 

points,  as  well  as  the  root  location  with  no  feedback,  may  be  calculated. 
The  damping  increments  as  calculated  from  the  locus  tue 

From  K*  =  O.OM*  second  to  K*  «  0.030  second 

V  B 

=  +  U.l  per  second 

From  K*  =  O.OUU  Second  to  K*  =  0.060  second 

e  e 

AMg  =  -  2.3  per  second 


Now  the  location  of  the  unaugmented  roots  of  the  vehicle  may  be  calculated 
by  finding  the  root  location  where  A Mg  *  ♦  8.15  per  second.  This  calcu¬ 
lation  yields  for  the  characteristic  roots  of  the  unaugmer.ted  motion  (  ©  ) 


s  =  1.00  ±  l.l*5i  per  second 


The  transient  motion  corresponding  to  this  calculated  result  agrees 
closely  with  the  time  history  of  the  measured  model  motion  shown  in 
Figure  18.  Ihe  0-wf  motion  is  shown  in  Figure  2h.  This  motion  was  not 
analyzed  due  to  the  highly  unstable  character  of  the  motion. 

Now  we  proceed  to  consider  the  three-degree-of-freedom  motion,  using  the 
data  with  feedback  shown  in  Figure  28.  Data  at  other  levels  of  feedback 
are  shown  in  Figures  26  and  27.  It  is  assumed,  as  indicated  by  the  static 
isolated  duct  data  of  Reference  1,  that  is  negligible.  In  this  case, 
the  complete  characteristic  equations  (as  obtained  from  equations  (15), 
where  Xy.  =  0,  z =  0,  xpg  =  0,  and  kg  =  0)  are 


‘  Vr  *  it  '  Vr  -  Vo  9 


Mu“r  -  -  Vf  +  9  -  (Me  *  Vof>  9  -  Vof  9 


0 

0 

0 


.  (3M 
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The  characteristic  equation  can  be  arranged  in  the  following  form: 


where  a  ia  the  characteristic  equation  for  the  two-degree- of- freedom 

6  ,u 

motion,  and  (g  )  is  the  characteristic  equation  of  the  two-degree-of- 

9  *u  o 

freedom  motion  with  ■  0.  The  coefficients  of  these  two  polynomials  may 
be  calculated.  is  known  from  the  one -degree -of- freedom  analysis,  and 
so  we  may  find  the  two  unknowns  2^  and  Zy  as  follows  by  placing  this 
equation  in  root  locus  form,  considering  and  Zy  as  unknown. 


The  poles  of  this  expression  arc  known,  having  been  determined  from  the 
two-degree -of- freedom  analysis.  Sot#  that  this  case  has  a  feedback  gain 

■  0.030  second 

Thus,  the  characteristic  equation  that  determines  the  poles  is 

Vu  "  (*  *  °’936  Xu)(*“  *  (fte  4  HftU°f)  1  “  4  0,936  8  (37) 

Substituting  the  previously  obtained  values  of  the  derivatives,  this 
expression  : -educes  to 


a  ,  -  s*  ♦  5.67  s*  -  2.91  «  ♦  11.83  (  38) 

v  »u 

The  roots  of  this  polynomial  are 

•1  •  -  6.29  per  second 

■  ♦  0.31  ±  1.251  per  second 

and  are  the  poles,  along  with  s  ■  0,  of  the  root  locus  expression, 
equation  (36),  and  are  shown  on  Figure  29  ( © ). 
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Part  of  the  numerator  of  equation  (36)  is  calculated  with  My  »  0  and  is 
equal  to 


(A*  „)  -  s3  ♦  5.67  ♦  3.0U  ♦  15. U 

W  »U  o 

The  roots  of  this  polynomial  are 

st  ■  -  5.18  per  second 

*  *  O.O65  *  1.731  per  second 

\ 

and  are  also  shewn  on  Figure  29  (O')* 

First  **  sketch  the  root  locus  for  the  factor  in  brackets  in  the  numerator, 
expressing  it  in  root  locus  form  as 


u) 

2»u  0 


(39) 


Nov,  this  root  locus  is  drawn  on  Figure  29  (•*•),  and  the  location  of  roots 
(taken  in  conjunction  with  the  complete  root  locus  equation),  that 
satisfies  the  l80°-angle  condition  at  the  three -degree -of -freedom  experi¬ 
mental  roots  (  Q  )  is  determined.  The  e xpe r ism nt ally  determined  charac¬ 
teristic  roots  for  the  three -degree -of -freedom  motion  are. 


s,  ■  ♦  0.25  *  1.321  per  second 

The  location  of  the  roots  of  equation  (39),  which  are  the  zeros  of  the 
root  locus  equation  (36),  was  found  to  be  at 

S|  ■  -  U.56  per  second 

s|()  ■  -  0.32  a  0.761  per  second 

and  are  shown  on  Figure  29  as  (  Q  ). 

Then  the  gain  at  this  point  is  calculated  to  yield  the  value 

Z,,g 

My  — -  ■  13.1  per  second  cubed  (Uo) 


Now  we  may  sketch  the  complete  root  locus  based  o»~.  equation  (36)  to 
determine  Zy  with  poles  (  O  )  •*»<!  zeros  (  Q  ).  The  locus  is  drawn  with 

Zy  varying  and  the  ratio  Zy/Zy  constant  as  shown  on  ®*igure  29  (  — ). 
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The  gain  calculation  yields 

2^  =  -  0.137  per  second 

Then,  cnowing  M^,  the  value  of  ^  is  calculated  from  equation  (^0)  as 

Zy  =  -  0.19  per  second 

Hiis  completes  the  analysis  of  the  70°  case  and  this  numerical  example. 

The  analysis  of  the  other  cases  follows  a  similar  procedure  as  described 

in  the  text. 
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Figure  17.  Self-Excited  Transient  Response.  One  Degree  of  Freedom,  9. 
No  Stability  Augmentation,  i^  =  70°,  0,7B(I  ■  25.2°. 
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Figure  21.  Self-Excited  Transient  Responses.  IVro  Degrees  of  Freedom, 
0-uf.  Kg  =  0.030  sec.  id  =  70°,  0<7BR  =  26.2°,  rpm  =  6780 
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Figure  22.  Self -Excited  Transient  Responses.  Two  Degrees  of  Freedom, 
0-Uf.  K*  =  0.044  sec.  id  =  70°,  p#78„  =  26.2°,  rpm  =  6780. 
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Figure  23. 


Self-Excited  Transient  Responses.  Two  Degrees  of  Freedom, 
6-Uf.  Kg  -  0.060  sec.  id  -  70°,  e>7#,  -  26.2°,  rpm  -  6780. 
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F_gure  25.  Self -Excited  Transient  Responses.  Three  Degrees  of  Freedom, 

Ic  Stability  Augmentation.  id  •  70°.  0#78B  •  26.2°, 

rptn  =  678c. 
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Figure  26.  Self-Excited  Transient  Responses.  Three  Degrees  of  Freedom, 

e-Uf-v-.  »  0.021  sec.  id  •  70°,  B#7#t  ■  26.2°,  rpm  »  6700. 

j2 


Figure  29.  Continued 


APPENDIX  III 

CONVERSION  TO  FULL  SCALE 


The  results  of  the  model  experiments  may  be  converted  to  cormsp  .  1  to  a 
variety  of  full-scale  vehicles  of  similar  geometry  to  the  model.  ±t  is 
convenient  to  consider  the  scaling  of  the  data  in  two  steps. 

First,  the  size  of  the  model  is  accounted  for  by  using  the  dynamic  model 
scale  factors  given  in  Table  II.  The  full-scale  aircraft  of  interest  here 
has  a  linear  scale  factor  of  0.145;  however,  other  scale  factors  may  be 
selected  to  correspond  to  ether  geometrically  similar  aircraft  of  desired 
size.  This  scaling  will  imply  a  certain  gross  weight  for  the  full-scale 
vehicle. 

Second,  the  results  may  be  interpreted  at  other  gross  weights  by  varying 
certain  of  the  parameters  involved,  maintaining  the  lift  coefficient  (or 
equivalently,  the  propeller  thrust  coefficient  based  on  forward  speed'' 
constant.  As  the  gross  weight  is  varied,  either  the  forward  speed  or  the 
ambient  air  density  can  be  varied  to  preserve  the  equilibrium  lift 
coefficient. 

These  two  interpretations,  and  the  appropriate  factors  to  use  for  gross 
weight  variation,  are  given  in  Table  V.  We  consider  here  only  the  effects 
cf  changes  in  gross  weight;  the  size  considerations  have  been  taken  into 
account. 

VELOCITY  -  r.ROSS  WEIGHT  CORRESPONDENCE 


Maintaining  the  equilibrium  lift  coefficient  of  the  vehicle  at  two  differ¬ 
ent  gross  weights,  at  the  same  altitude,  yields  the  following  relationship 
between  flight  velocity  and  gross  weight: 

4  \ 


Defining  a  weight  ratio  scale  factor  as 


Aw 


the  velocity  la  scaled  as 


n 


The  advance  ratio  must  also  be  maintained  constant,  and  so  this  scaling 
results  in  a  different  rpm;  i.e., 


In  the  case  of  the  experiments  conducted,  0B  will  not  correspond  to  the 
full-scale  rotational  speed  of  the  ducted  propeller,  since  the  model  rpm 
was  selected  on  the  basis  of  a  proper  value  of  .  Scale  factors  for 
conversion  of  the  data  in  this  fashion  are  given  in  Table  V. 

It  is  possible  to  make  an  assumption  that  will  make  a  wider  interpretation 
of  the  data  possible.  This  assumption  is  similar  to  the  use  of  the  pro¬ 
peller  thrust  coefficient  to  characterize  data  on  tilt-wing  aerodynamics. 
That  is,  if  we  assume  that  the  aerodynamic  stability  characteristics  of 
the  vehicle  depend  primarily  on  the  ratio  of  forward  speed  to  duct  exit 
velocity  and  not  on  the  particular  combination  of  blade  angle  and  advance 
ratio  used  to  produce  this  velocity,  it  may  be  assumed  that  blade  angle 
and  rpm  are  interchangeable  and  the  scaled  data  may  be  applied  to  other 
rotational  speeds. 

The  validity  of  this  assumption  has  not  been  checked.  Hovering  flight 
data  indicate  that  there  are  differences  in  the  dynamics  depending  upon 
the  combination  of  blade  angle  and  rpm  used  to  produce  a  given  thrust,  so 
this  approximation  should  be  applied  with  care  if  used  at  all. 

AIR  DENSITY  -  GROSS  WEIGHT  CORRESPONDENCE 

Alternately,  the  lift  coefficient  may  be  maintained  constant  by  varying 
ambient  air  density  in  proportion  to  gross  weight: 

K  _  ^ 

*  "  Po 

Then  the  data  may  be  interpreted  on  this  basis,  where  the  aerodynamic 
forces  will  vary  by  the  scale  factor  \r  and  the  reduced  gross  weight  will 

be  equivalent  to  flight  at  a  different  altitude  given  by 


In  this  cue,  note  that 


The  scale  factors  for  conversion  by  this  method  are  also  given  in  Table  V. 
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In  this  case,  it  may  be  observed  that  there  will  be  no  change  in  the  dy¬ 
namic  stability  characteristics  of  the  aircraft.  This  result  indicates 
that  in  many  cases  it  is  desirable,  for  comparison  with  flight  test,  to 
test  a  model  that  is  overweight  on  the  basis  of  the  dynamic  scaling  law, 
since  the  flight  test  experiments  will  always  be  conducted  at  altitudes 
above  sea  level.  This  correspondence  for  the  experiments  conducted  here 
is  shown  in  Figure  lU. 
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discussed. in  this  report. 


The  measured  time  histories  indicated  that  the  data  could  be  analyzed  on  the  basis  of 
linearized  small  perturbation  equations.  Root  locus  techniques  were  used  to  analyze 

the  data. 

The  full-scale  derivatives  determined  from  the  analysis  that  correspond  to  a  vehicle 
very  similar  to  the  Bell  X-22A  are  presented. 
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the  lowest  incidence  investigated. 
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